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Abstract: Ecological biogeography of small mammal assemblages in the central Andes of Argentina. 

The Andes constitutes an extensive mountain range and a true laboratory for the study of the evolution of South 
American biota. In recognition to José Yepes, one of the pioneers in characterizing the distribution of mammals of 
Argentina we seek in this contribution to synthesize, from ecological biogeography, the knowledge we have about 
the diversity patterns of small mammals of the Central Arid Andes (ACA). We characterize regional and local 
patterns of rodent richness and endemism as a function of latitude, altitude and area. At regional scale, we con- 
structed a species database using species range maps, meanwhile at the local scale, we conducted small mammal 
surveys using standardized techniques along 4 elevational transeets between 32 Q S and 35 Q S. Among the results, 
we highlight, at regional level: a) a high number of endemisms (> 50%); b) richness increase with elevation, and 
c) endemic species - area relationship. Among the results at local level we stress: d) richness increases with el¬ 
evation, mean precipitation and topographic heterogeneity, and e) greater abundance at intermediate elevations. 
Biogeographic-ecological research allows us to depict patterns of composition and distribution of the central 
Andean biodiversity, and assess some of its causal factors. However, the long term conservation of the Andes re- 
quires not only solid biological / ecological knowledge but also conservation policies that intégrate the biological, 
social and cultural dimensions of the Andean ecosystem. This approach seems the most appropriate alternative 
in the face of the different impaets induced by man and the climate change in the región. 

Key words: Rodents, Andes, environmental gradient, species diversity. 

Resumen: Los Andes constituyen una de las cadenas montañosas más extensas del planeta y un verdadero labo¬ 
ratorio para el estudio de la evolución de la biota sudamericana. En homenaje a José Yepes, uno de los pioneros 
en caracterizar la distribución de mamíferos de Argentina, buscamos en esta contribución sintetizar, desde la 
biogeografía ecológica, el conocimiento que tenemos sobre los patrones de diversidad de pequeños mamíferos de 
los Andes Centrales Aridos (ACA) y los factores causales. Caracterizamos, a nivel regional y local, los patrones de 
riqueza de roedores y endemismos en función de latitud, altitud y área. A escala regional, generamos una base de 
datos a partir de rangos de distribución de las especies, mientras que a escala local se realizaron muéstreos en 4 
transectas altitudinales entre los 32 Q S y 35 Q S. Entre los resultados destacamos, a escala regional: a) un elevado 
número de endemismos (> 50 %); b) aumento de riqueza con la altura, y c) relación especies endémicas - área. A 
escala local destacamos: d) la riqueza incrementa con la altura, precipitación media y heterogeneidad topográfica, 
y e) mayor abundancia en alturas intermedias. 

Palabras clave: Roedores, Andes, gradientes ambientales, diversidad. 


INTRODUCCIÓN 

Sobre la base de los estudios de Oldfield 
Thomas y destacando el valor de las coleccio¬ 
nes científicas como base para diversos estudios 
taxonómicos y de faunas regionales, José Yepes 
es uno de los pioneros, junto a Angel Cabrera, 
de caracterizar la distribución geográfica de los 
mamíferos de Argentina (1930, 1937,1938, 1941; 
véase García, 2018 para una síntesis biográfica de 


Yepes) y el primero en sintetizar la zoogeografía 
de mamíferos de la provincia de Mendoza (1936). 
No obstante el limitado relevamiento regional, 
escaso número de ejemplares y conocimiento 
taxonómico de un reducido número de especies, 
Yepes caracteriza a los mamíferos de Mendoza a 
modo de “catálogo documentado”, según sus pa¬ 
labras, dentro de 5 divisiones o provincias zooló¬ 
gicas propuestas por Lahille (1898), una de las 
cuales es la “provincia andina”, o equivalente a 
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la región montañosa de los Andes centrales que 
tratamos en el presente trabajo. 

Los ecosistemas montañosos han desempe¬ 
ñado un papel importante en la evolución y dis¬ 
tribución de la biodiversidad (Lomolino, 2001; 
Kórner, 2000). La orogenia de la Cordillera de los 
Andes tuvo un impacto substancial en el desarro¬ 
llo y evolución biótica de la Región Neotropical, 
imponiendo barreras físicas y limitantes altitudi- 
nales a la distribución de numerosos taxa y a la 
migración de distintos linajes (Patterson et al., 
2012 ). 

La influencia climática de la orogenia Andina 
ha determinado cambios en el paisaje, originando 
un mosaico de nuevos ambientes y ecosistemas 
áridos (ej. Patagonia y Puna), con una amplia 
variedad de relieves topográficos, climas y hábi¬ 
tats. Esto contribuyó al desarrollo de una gran 
diversidad biótica en la cordillera de los Andes 
y mantenimiento de numerosos endemismos 
(Hershkovitz, 1969; Patterson & Pascual, 1972; 
Reig, 1981, 1986; Marquet, 1994). El levanta¬ 
miento de los Andes, promovió el aislamiento de 
organismos de tierras bajas a ambos lados y en 
los sectores más altos de la cordillera (Marshall, 
1979; Palma, et al., 2005; Brumfield & Edwards, 
2007), y actuó como corredor norte - sur y sur- 
norte durante el gran intercambio biótico-ameri- 
cano (Gabi -Marshall, 1979). 

Los patrones de diversidad a escala global y 
regional a menudo se asocian con gradientes de 
latitud y elevación (Simpson, 1964; Terborgh, 
1977; Rahbek, 1995, 1997; Rosenzweig, 1995; 
Brown, 2001; Lomolino, 2001). Así, la disminu¬ 
ción de la diversidad a medida que aumenta la 
latitud es uno de los patrones más generalizados. 
Este ha sido registrado en numerosos taxa, ac¬ 
tuales y extintos, y a diferentes escalas, desde lo¬ 
cal a global (Willig et al., 2003; Willig & Presley, 
2013). La latitud puede ser considerada un subs¬ 
tituto de varios factores, entre los que se destaca 
el área, el clima, la productividad, la heterogenei¬ 
dad topográfica y la altura (Willig & Bloch 2006; 
Badgley & Fox 2000; Currie et al., 2004). En el 
caso de los mamíferos, la forma y fuerza de este 
patrón varía según el grupo de mamíferos estu¬ 
diados (Mares & Ojeda, 1982; Ruggiero, 1994, 
Ruggiero et al., 1998). 

Los gradientes altitudinales de diversidad 
constituyen uno de los patrones biogeográfi- 
cos mas conocidos y estudiados, remontándo¬ 
se a los trabajos de Darwin (1839, 1859) en los 
Andes centrales de Chile; Wallace (1876, 1878) 
en Indonesia; Merriam (1890); Whittaker (1960); 
Whittaker & Niering (1965) y Brown (1971) para 


Norte América, entre otros. Actualmente se co¬ 
nocen 5 patrones generales sobre la relación de 
la diversidad de especies con la elevación: 1) el 
número de especies disminuye monotónicamen- 
te con la elevación; 2) riqueza alta en la porción 
inferior del gradiente que se mantiene en for¬ 
ma de meseta, seguida por una disminución de 
la riqueza; 3) riqueza alta en elevaciones bajas 
con un máximo de diversidad encontrado en al¬ 
turas intermedias; 4) pico de riqueza de especies 
en altura intermedia, es decir un pico unimodal 
de diversidad. 5) la riqueza de especies aumen¬ 
ta con la elevación (Martin, 1958; Wake et al., 
1992; Grytnes et al., 2007; Novillo & Ojeda, 2012, 
Novillo & Ojeda, 2014). Algunas de las posibles 
explicaciones para estos patrones, establecen que 
los sitios de mayor diversidad se relacionan con 
regiones de condiciones ambientales óptimas, o 
bien pueden registrarse en zonas donde las espe¬ 
cies se superponen, o en lugares donde comuni¬ 
dades de vegetación distintas se encuentran en 
estrecha proximidad (Lomolino, 2001). 

Los factores causales de los patrones de diver¬ 
sidad pueden agruparse en 4 categorías: clima, 
espacio, historia evolutiva y procesos ecológicos 
(Pianka, 1966; Gastón, 2000; McCain, 2007). 
Las hipótesis climáticas se basan en la variación 
de variables abióticas como temperatura, preci¬ 
pitación y humedad, entre otras. Las hipótesis 
espaciales incluyen la relación especie-área clási¬ 
ca (SAR). Esta expone un aumento de riqueza a 
medida que aumenta el área (Rosenzweig, 1995). 
Y puede explicarse mediante dos mecanismos: 1) 
la probabilidad de extinción es menor en áreas 
mayores (MacArthur & Wilson, 1967), y 2) áreas 
mayores contienen mayor diversidad de ambien¬ 
tes, favoreciendo la coexistencia de un mayor nú¬ 
mero de especies (Simberloff, 1976). La historia 
evolutiva, se refleja en la variación de las tasas 
de especiación, extinción y edad de los ciados. Y 
las hipótesis ecológicas incluyen, los efectos eco- 
tonales, la competencia, la heterogeneidad y la 
complejidad del hábitat, entre otros. También es 
necesario remarcar que los patrones de diversi¬ 
dad y sus factores causales varían en función del 
tamaño de rango de distribución de las especies. 
En general, especies con rangos geográficos pe¬ 
queños tienden a mostrar un patrón de riqueza 
diferente al de las especies con distribuciones 
más extendidas (Jetz et al., 2004). 

Los estudios sobre gradientes de diversidad 
a lo largo de los Andes se han centrado princi¬ 
palmente en ambientes tropicales y húmedos 
(Patterson et al., 1998; Brehm & Fiedler, 2003; 
Kattan & Franco, 2004; Herzo get al., 2005, Ferro 
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& Barquez, 2009; Ferro, 2013, etc), destacando la 
relación especie - energía (en un sentido amplio, 
teniendo en cuenta temperatura y productivi¬ 
dad- Wright, 1983), y el equilibrio hídrico como 
principales determinantes de los gradientes de 
diversidad (Bini et al., 2004). Estudios similares 
sobre la porción árida de los Andes son relativa¬ 
mente escasos y la mayoría se refieren al desierto 
de altura de la Puna (Pearson & Pearson-Ralph, 
1978; Marquet, 1994). 

En este trabajo caracterizamos, por un lado, 
los patrones de diversidad de especies y endemis- 
mos de roedores en función de latitud, altitud y 
área en los Andes áridos centrales, tanto a nivel 
regional, entre los 8 o y 42° S (Young et al., 2007; 
Sempere et al., 2008) como a nivel local (Provincia 
de Mendoza). Por otro lado, analizamos los facto¬ 
res topográficos, climáticos y espaciales involu¬ 
crados en el modelado de estos patrones. 

METODOLOGÍA 

Para el análisis de los patrones de distribu¬ 
ción de las especies de roedores andinos, a escala 
regional, se utilizaron 2 aproximaciones: 1) se 
analizaron los patrones de distribución para to¬ 
das las especies presentes en la región de estu¬ 
dio propuesta (ensamble total); 2) se utilizaron 
solo las especies endémicas. Para la escala local 
se consideró el pool de especies registradas en 
muéstreos de campo. 

Escala regional 

Area de estudio: El área de estudio propuesta 
como escala regional, se encuentra en el segmen¬ 
to central de los Andes (desde los 8°S hasta los 
42° S) (Fig. 1A) y se denomina Andes Centrales 
Aridos (ACA). Presenta una extensión altitudinal 
desde los 800 hasta los 6000 msnm. Para gene¬ 
rar una base de datos de ocurrencia de pequeños 
mamíferos pertenecientes a los ACA, se utiliza¬ 
ron las áreas de distribución confeccionadas en 
la base de datos Natureserver (Patterson et al., 
2007), las cuales fueron ajustadas mediante bi¬ 
bliografía actualizada. Se utilizó literatura cien¬ 
tífica general y particular de la región de estu¬ 
dio (Pearson, 1982; Redford & Eisenberg, 1992; 
Marquet et al., 1998; Díaz & Barquez, 2002; 
Wilson & Reeder, 2005). Los rangos de distribu¬ 
ción se utilizaron para confeccionar una matriz 
de presencia y ausencia por bandas latitudinales 
y altitudinales. El rango latitudinal se dividió en 
bandas de 2 grados y el altitudinal en bandas de 
200 m (Hurlbert & Jetz, 2007). Se cuantificó el 
número de especies presentes en cada banda para 


examinar la relación entre el número de especies 
y la latitud-altitud. Se documentaron patrones 
de riqueza de especies a lo largo del gradiente 
usando modelos lineales generalizados. Cada mo¬ 
delo lineal generalizado se comparó utilizando el 
criterio de información de Akaike (AIC). 

El área por banda altitudinal se calculó a par¬ 
tir de un modelo de elevación (Jarvis et al., 2008). 
Se realizó un proceso de reclasificado de este 
DEM en bandas de 200 m de altura y para cada 
banda se calculó el área mediante la herramienta 
“zonal statistics” de ArcGis 9.1. Para calcular el 
efecto del área a lo largo de los ACA, se realiza¬ 
ron regresiones lineales log-log entre la riqueza 
y el tamaño del área (latitudinal y altitudinal). 
Calculamos los valores de riqueza corregidos por 
el tamaño del área utilizando el método curvilí¬ 
neo (c = S / A z , siendo c = estimación de la diver¬ 
sidad de área corregida, S = riqueza de especies, 
A = tamaño del área y z = factor constante. Se 
utilizó un valor z específico de taxón para peque¬ 
ños mamíferos no voladores (z = 0,22) calculado 
por McCain (2007). 

Escala local 

El estudio se realizó en los Andes centrales de 
Argentina, ubicados entre los 30° S y 36° S en la 
provincia de Mendoza (Fig. IB). El área localiza¬ 
da en las ecoregiones de los Altos Andes y Monte, 
se caracteriza por una vegetación que varía 
desde matorrales, en las estribaciones andinas 
(1300 m sobre el nivel del mar), hasta arbustos 
enanos, gramíneas y hierbas, en elevaciones más 
altas (3300 m sobre el nivel del mar) (Martínez 
Carretero & Méndez, 1992; Mendez et al., 2006). 

Se establecieron 4 transectas de muestreo, 
ubicadas entre los 30° y 36° de latitud sur (TU, 
TT, TD, VH; Fig. Ib). Cada transecta presenta un 
rango altitudinal desde los 1300 a 3300 msnm. 
En cada una de ellas se seleccionaron 5 sitios de 
muestreo a distintas alturas (1300, 1800, 2300, 
2800, 3300 msnm). En cada sitio se colocaron 
150 trampas tipo Sherman de captura viva, abar¬ 
cando la mayor cantidad de ambientes posibles. 
Cada sesión de trampeo consistió en 3 noches por 
sitio. Las trampas fueron cebadas con pasta de 
maní y avena arrollada. El estudio se desarrolló 
durante dos años consecutivos, relevando todas 
las transectas en dos oportunidades. El esfuerzo 
de muestreo total fue de 18000 noches trampa. 

Para la caracterización de la estructura de 
los ensambles (riqueza, abundancia y equitati- 
vidad) utilizamos curvas de rango abundancia 
(Whittaker, 1972). En estas curvas las especies 
se grafican en secuencias/rangos desde las más 
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Fig. 1. A) Mapa de la ubicación geográfica de los Andes Centrales Áridos. B) Ubicación de la zona de estudio en 
los Andes centrales de Argentina, Provincia de Mendoza. Se registra el gradiente en elevación y las 4 transectas 
altitudinales. TU: transectaUspallata; TT: transectaTunuyán; TD: Transecta Laguna del Diamante; VH: tran- 
secta Valle Hermoso. 


abundantes a las menos abundantes a lo largo 
del eje X. A su vez sobre este eje se puede ob¬ 
servar el número de especies presentes en el 
ensamble, cuanto más extenso, mayor riqueza 
presenta el ensamble. Estas curvas, también nos 
dan información sobre que modelo de abundan¬ 
cia describe mejor nuestros datos. Las curvas 
muy empinadas representan ensambles con alta 
dominancia (distribuciones geométricas o series 
logarítmicas), mientras que curvas con menor 
pendiente (distribución “log normal” o en “vara 
quebrada”) dan cuenta de una mayor equitativi- 
dad (Magurran, 2004). 

Para cuantificar la variación de riqueza de 
pequeños mamíferos a lo largo del gradiente 
altitudinal se utilizo el número de especies por 
banda altitudinal. La relación entre la riqueza y 
la elevación se evaluó utilizando modelos lineales 
generalizados y mediante análisis polinómico. El 
modelo de mejor ajuste se seleccionó utilizando 
el criterio de información de Akaike (AIC) y, para 
los modelos que presentaron diferencias AIC me¬ 
nores de 2, se consideró el criterio R 2 . 

El efecto del clima y la topografía sobre la 
diversidad de pequeños mamíferos se evaluó 
utilizando modelos lineares generalizados me¬ 
diante selección por diferencias de AIC y criterio 
R 2 (Crawley, 1993), utilizando el software SAM 
(“Spatial Analysis in Macroecology”, Rangel et 
al., 2010). 

Para separar los efectos conjuntos e indepen¬ 
dientes de clima, topografía y espacio, se llevo a 


cabo un análisis de partición de la varianza. Este 
análisis se llevó a cabo mediante el uso de regre¬ 
siones parciales (Legendre & Legendre, 1998, 
Lobo et al., 2002, Hortal et al., 2008). Utilizamos 
la ecuación polinomial de tercer grado (análi¬ 
sis de superficies de tendencias - Legendre & 
Legendre, 1998) como un proxy para los descrip¬ 
tores espaciales. Estos análisis se realizaron en el 
software R, versión 2.15.2 (R Development Core 
Team 2012). 

RESULTADOS 

Escala regional 

A lo largo de los ACA se registra un total de 
sesenta y siete especies de pequeños mamíferos 
no voladores, pertenecientes a 7 familias y 26 
géneros (Apéndice 1). Treinta y tres (49%) son 
endémicas del área de estudio propuesta. 

La riqueza de especies a lo largo del gradiente 
latitudinal mostró un patrón bimodal (Fig. 2A), 
un pico de riqueza se registró a 18° S con 34 es¬ 
pecies y un segundo pico de 25 especies ocurrió 
a los 36° S. Se encontró el mismo patrón bimo¬ 
dal para especies endémicas, mostrando menos 
especies en altas latitudes. El patrón altitudinal 
de riqueza mostró un pico entre los 3000 y 3500 
msnm. Las especies totales presentaron una for¬ 
ma de joroba no simétrica, con una riqueza de 
especies decreciente en altas elevaciones. Por 
otra parte, las especies endémicas presentaron 
una forma simétrica con una disminución en la 
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Fig. 2. A) patrón latitudinal de variación de la riqueza de especies. B) patrón altitudinal de variación de la ri¬ 
queza de especies. C) Distribución de la riqueza y el área a lo largo del gradiente latitudinal.D) Distribución 
de la riqueza y el área a lo largo del gradiente altitudinal. E) Relación Log-log entre riqueza de especies y área 
por banda latitudinal. F) Relación Log-log entre riqueza de especies y área por banda altitudinal. G) Relación 
entre la riqueza corregida por área (c) y la latitud. H) Relación entre la riqueza corregida por área (c) y la al¬ 
titud. Círculos negros: todas las especies. Cuadrados blancos: especies endémicas. Triángulos negros: área. 
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Fig. 3. Curvas rango abundancia para cada transecta altitudinal. A) Transecta Las Leñas-Valle Hermoso (VH); B) 
Transecta Laguna del Diamante (TD); C) TransectaTunuyán (TT); D) TransectaUspallata (TU). 1) E. morgani, 
2) P xanthopygus, 3) Phyllotissp, 4) A.olivacea, 5) E. chinchilloides, 6) L. micropus, 7) A. andina , 8) C. macronyx, 
9) Akodon spegazzinii, 10) E. mordax, 11) G. griseoflavus; 12) E. moreni; 13) C. musculinus; 14) A. dolores. Log 
(pi)= species relative abundace. 


riqueza de especies hacia ambos extremos del 
gradiente (Fig. 2B). 

El tamaño del área varía a lo largo de ambos 
gradientes. A los 30° S, se detectó una fuerte 
disminución en el tamaño de área, que coincidió 
con una disminución en la riqueza de especies 
para ambos conjuntos de datos (especies totales 
y endémicas, Fig. 2C). El tamaño de área para 
las bandas latitudinales representó el 25% de la 
variación de la riqueza para las especies totales 
(. R 2 = 0,25; p = 0,04) y el 70% para las especies 
endémicas (R 2 = 0,70; p <0,001) (Fig. 2E). 

Cuando se utilizo la estimación de riqueza 
corregida por el área (c = SA Z ) la relación con la 
latitud fue positiva para las especies totales (R 2 
= 0,35; p = 0,014) y ninguna tendencia significa¬ 
tiva para las especies endémicas (R 2 = 0,0 6; p = 
0,3) (Fig. 2G). 

A lo largo del gradiente altitudinal el tamaño 
del área presentó sus valores máximos en alturas 
intermedias (3000-4500 m) (Fig. 2D). La relación 
entre la riqueza y el área representó el 90% de 
la variación del gradiente altitudinal para el con¬ 
junto de especies totales (R 2 = 0,91; p <0,001) y 
el 20% para el conjunto endémico (R 2 = 0,20; p = 
0,01) (Fig. 2F). Cuando la riqueza de especies se 
ajustó por tamaño de área, no hubo relación sig¬ 
nificativa a lo largo del gradiente altitudinal para 
el grupo de especies total (R 2 = 0,06; p = 0,18). 


Mientras que para las especies endémicas la re¬ 
lación fue positiva y significativa, representando 
el 57% de la variación de riqueza con altitud (R 2 
= 0,57; p <0,001) (Fig. 2H). Los polinomios de 
tercer orden fueron considerados los mejores mo¬ 
delos para explicar los gradientes latitudinales y 
de elevación (Tabla 1). 

Escala local 

A escala local cuantificamos 9 géneros y 14 
especies de roedores sigmodontinos (. Abrothrix 
olivácea, Abrothrix andina, Akodon spegazzinii, 
Akodon dolores, Phyllotis xanthopygus, 
Calomys musculinus, Chelemys macronyx, 
Euneomys chinchilloides, Euneomys mordax, 
Graomys griseoflavus, Loxodontomys micropus, 
Eligmodontia moreni, Eligmodontia morgani 
y Phyllotis sp ). Se capturó un total de 1833 
individuos con un éxito de captura total del 
10%; el mayor éxito de captura (16% a 29%) se 
registró en alturas intermedias (entre 1800 y 
2300 msnm). 

Cada transecta altitudinal presentó 
diferentes patrones comunitarios, en cuanto a la 
composición de especies, riqueza y equitatividad 
de los ensambles. A lo largo de la transecta Las 
Leñas - Valle Hermoso (VH) (35° S) se registraron 
10 especies de pequeños mamíferos (Fig. 3A). 
Las bandas altitudinales inferiores (1300 - 
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Tabla 1. Relación entre la riqueza corregida por área (c) y la latitud y altitud. 



Modelo 

AIC 

dAIC 

Peso 


c-Ts 

c Es 

c-Ts 

c Es 

c-Ts 

c Es 

Patrón latitudinal 

Linear 

46.8 

46.4 

4.1 

22.5 

0.06 

<0.001 


Cuadrático 

43.1 

33.6 

0.4 

9.7 

0.41 

0.007 


Cubico 

2.7 

23.9 

0 

0 

0.51 

0.99 

Patrón altitudinal 

Linear 

35.3 

22.4 

19.2 

18.6 

<0.001 

<0.001 


Cuadrático 

27.7 

12.5 

11.6 

8.8 

0.003 

0.01 


Cubico 

16.1 

3.7 

0 

0 

0.99 

0.98 


Comparaciones entre modelos lineales, cuadráticos y cúbicos. La negrita indica el mejor modelo. 
c-Es: especies endémicas; c-Ts: especie total;AIC:criterio de información de Akaike,; dAIC: delta AIC; 
desviación; peso: pesos Akaike. 


1800 msnm) mostraron menos especies que las 
más altas (2300 - 2800 - 3300 msnm). Por otro 
lado, las bandas más bajas presentaron mayor 
equitatividad que las más altas. La dominancia 
de las especies también varió con la altitud; E. 
morgani se registró como especie dominante a 
1300 (msnm), mientras que A. olivácea dominó 
entre 1800 y 2300 (msnm); C. macronyx lo es a 
2800 (msnm) y E. chinchilloides a 3300 (msnm). 

A lo largo de la transecta de Laguna del 
Diamante (TD) se registraron 10 especies que 
pertenecen a 7 géneros (Fig. 3B). La mayor 
riqueza correspondió con los sitios más elevados 
2800 - 3300 (msnm). La dominancia de las 
especies también varió con la altitud; E. morgani 
fue dominante en sitios entre 1300 - 1800 
(msnm), mientras que A. olivácea fue dominante 
en los sitios entre 2300 - 3300 (msnm). 

La transecta de Tunuyán (TT) exhibió 
9 especies de 7 géneros. La mayor riqueza y 
equitatividad de especies se registró entre los 
2800 y 3300 msnm (Fig. 3C). En esta transecta, 
las especies dominantes difirieron para cada 
banda altitudinal (E. moreni a 1300 msnm, P. 
xanthopygus a 1800 y 2300 msnm, A. spegazzinii 
a 2800 msnm y A. andina a 3300 msnm). 

En la transecta Uspallata (TU) se registraron 
8 especies y 7 géneros. Los sitios intermedios 
registraron valores de riqueza más bajos, 
mientras que los 1300 (msnm) y los 3300 (msnm) 
registraron mayor riqueza y equitatividad (Fig. 
3D). Los sitios más bajos fueron dominados por 
C. musculinus 1300 (msnm) y E. moreni 1800 
(msnm). P xanthophygus fue la especie dominante 
en las bandas intermedias 2300y 2800 (msnm); y la 
banda altitudinal más alta 3300 (msnm) presentó 
a E. chinchilloides como especie dominante. 


La riqueza de especies de pequeños mamíferos 
mostró una relación monotónica positiva con la 
elevación (29,19% de varianza explicada) (Fig. 
4A). Mientras que el patrón de abundancia 
evidenció una relación unimodal, con la mayor 
abundancia en alturas intermedias (55,95% de 
varianza explicada) (Fig. 4B). 

Varios predictores explicaron el patrón de 
elevación observado para la riqueza de especies 
(Tabla 2). La riqueza de especies se asoció 
positivamente con la precipitación media anual 
y la heterogeneidad topográfica, y se asoció 
negativamente con temperatura mínima, 
disminuyendo a medida que la temperatura 
aumentaba (Fig. 5A). 

La partición de la varianza muestra que 
las variables seleccionadas explican un 84% 
de la variación de la riqueza. El efecto puro 
de cada subconjunto de predictores (clima y 
topografía) fue similar y moderado (8,4% y 13%, 
respectivamente). Mientras que el efecto del 
espacio representó el 7,4% de la variación de la 
riqueza. El efecto conjunto del clima y topografía 
explican un 35% de la variación de la riqueza, 
mientras que el efecto puro de la topografía 
explicó un 13%. El efecto combinado del clima y 
el espacio representó el 12% de la variación de la 
riqueza total (Fig. 5B). 

DISCUSIÓN 

Más del 50% de los pequeños mamíferos de 
los Andes Aridos Centrales de Sudamérica son 
endémicos. Este alto grado de endemismo apoya 
observaciones anteriores sobre el rol determinante 
de la formación de ambientes áridos a grandes 
alturas producidos por el levantamiento de 
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Tabla 2. Efectos climáticos y topográficos sobre la riqueza de especies. Selección de modelos basados 
en diferencias de AIC y en criterio R 2 . 


Modelo climático: riqueza~AMT+MTCM+T(sd)+TAR*TMDR+AP+P(cv)+PET+AET+NDVI 






Delta 


Variables 

TzVar 

R 2 

AICc 

AICc 

AICcwi 

AP 

1 

0.586 

53.18 

0 

0.105 

MTCM + AP 

2 

0.632 

53.992 

0.812 

0.07 

AMT + AP 

2 

0.623 

54.459 

1.278 

0.055 

AP + PET 

2 

0.615 

54.902 

1.722 

0.044 

AP + P(cv) 

2 

0.596 

55.872 

2.692 

0.027 

Modelo topográfico: riqueza—Slope(sd)+DEM(rg)+NDVI (sd) 








Delta 


Variables 

nVar 

R 2 

AICc 

AICc 

AICcwi 

Slope(sd) + DEM(rg) + DEM(rg)2 

3 

0.542 

61.956 

0 

0.292 

DEM(rg) + Slope(sd)2 + DEM(rg)2 

3 

0.54 

62.05 

0.094 

0.279 

DEM(rg) + DEM(rg)2 

2 

0.446 

62.175 

0.22 

0.262 

DEM(rg) + Slope(sd) 

2 

0.343 

65.569 

3.613 

0.048 


PMA: precipitación anual; TMMF: temperatura mínima del mes mas frió; PET: evapotranspiración 
potencial; P(cv): coeficiente de variación de la precipitación. Pendiente (sd): desvío estándar de la 
pendiente; Dem (rg): Rango de elevación.La negrita indica el mejor modelo. 



Fig. 4. A) Relación entre riqueza de especies a escala 
local y altitud. B) Relación entre la abundancia relativa 
y la altitud. 


los Andes, en los procesos cladogénicos de los 
pequeños mamíferos Sudamericanos (Reig, 
1981; Pearson, 1982; Marquet, 1994). La 
singularidad de la biota andina se atribuye a las 
elevadas tasas de especiación promovidas por 
condiciones ambientales particulares: ambientes 
topográficamente complejos, bandas altitudinales 
homeotérmicas angostas o pequeñas, y una 
distribución discontinua o en parches del hábitat 
(Rahbek & Graves, 2000). 

La forma bimodal del patrón latitudinal 
de riqueza, con sus dos picos a los 18° S y 34 - 
35° S, y la marcada disminución de la riqueza 
a los 30° S, constituyen un área de divergencia 
biogeográfica. Una posible explicación de la 
disminución en la riqueza a los 30° S, está dada 
por la reducción de área Andina desde los 20° 
S en dirección sur, y la presencia de los picos 
montañosos más elevados de la cordillera, que 
pudieron actuar como barreras efectivas para 
varios taxa (Lanzone, 2009). Esta disminución 
en el área se ve acompañada por una marcada 
disminución en la elevación de la cordillera hacia 
latitudes mayores a los 33° S, favoreciendo de este 
modo la colonización de especies de procedencia 
patagónica. 

La fuerte correlación positiva entre la riqueza 
de especies endémicas y la elevación de los Andes 
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Fig. 5: A) Determinantes climáticos y topográficos de la riqueza de pequeños mamíferos a lo largo del gradiente al- 
titudinal de los ACA (PMA: precipitación media anual; TMMF: temperatura mínima del mes más frió. B) Diagrama 
de partición de la varianza de la riqueza, en función del conjunto climático, el conjunto topográfico y el conjunto 
espacial de variables predictoras. T (topografía), C (clima), S (espacio), U (varianza no explicada por el modelo). 


centrales áridos, contradice el patrón general de 
disminución de riqueza a mayor altura (Stevens, 
1992; Rosenzweig, 1995; Gastón & Blackburn, 
2000). Este resultado concuerda con lo observa¬ 
do en un estudio a menor escala, en el bioma de 
la Puna (Marquet, 1994), y refuerza el rol de los 
Andes en el proceso cladogenético de los roedores 
Sudamericanos (Reig, 1981; Leite et al., 2014). 

A escala local el patrón altitudinal de rique¬ 
za es coincidente con el patrón reportado para 
las especies endémicas a escala regional. Este 
aumento de riqueza en función de la altitud, es 
poco frecuente, ya que en la mayoría de los estu¬ 
dios sobre gradientes de diversidad en roedores 
se menciona al patrón unimodal como generali¬ 
dad (Heaney, 2001; Rahbek, 2005; McCain, 2005, 
McCain & Grytnes, 2010). Los escasos trabajos 
que reportan una mayor riqueza en zonas altas, 
pertenecen a zonas áridas de los Andes (Pearson 
& Ralph, 1978; Marquet, 1994; Patterson et al., 
1989, 1998; Ruggiero & Kitzberger, 2004). La re¬ 
lación entre altura y riqueza de especies varía se¬ 
gún el tipo de ambiente, las características climá¬ 
ticas de los gradientes y la historia evolutiva de 
la especies que componen los ensambles. De este 
modo, la riqueza de especies es elevada en zonas 
de mayor altura en la cordillera árida, mientras 
que en general se registra un patrón unimodal en 
las zonas tropicales. 

La interacción entre clima y topografía se 
registra como mejor predictor para el patrón de 
riqueza de especies de pequeños mamíferos andi¬ 
nos. Se cree que la variación climática asociada 
con la heterogeneidad topográfica genera aisla¬ 


miento geográfico local, lo que puede promover 
la especiación a través de la alopatría (Ruggiero 
& Hawkins, 2008). La productividad medida 
como la precipitación media anual, es uno de 
los factores que explican el patrón de riqueza de 
roedores andinos. La hipótesis de productividad 
plantea que la energía (productividad primaria) 
limita la diversidad mediante la regulación de las 
cascadas tróficas. La productividad representa 
la cantidad de recursos disponible para los con¬ 
sumidores, por lo tanto, zonas de mayor produc¬ 
tividad, permitirán la coexistencia de un mayor 
número de especies (MacArthur, 1972; Terborgh, 
1977, Wright 1983, Turner et al., 1987). Además, 
cuando los factores climáticos se ven modulados 
o relacionados con cambios en la elevación y re¬ 
lieve, estos pueden generar condiciones climáti¬ 
cas locales (microclimas), pueden afectar el área 
del hábitat disponible (microhabitats) y la conec- 
tividad entre estos hábitats (Badgley, 2010), pro¬ 
moviendo una mayor disponibilidad de nichos y 
de recursos. 

La hipótesis de energía ambiente (temperatu¬ 
ra), predice que una mayor temperatura permiti¬ 
rá que los organismos inviertan menor cantidad 
de energía (ej. tiempo y esfuerzo) en termorre- 
gulación, favoreciendo de esta manera el pro¬ 
cesos de reproducción. Esto permitirá mayores 
tamaños poblacionales y menor riesgo de extin¬ 
ción, promoviendo una mayor riqueza de espe¬ 
cies (Turner et al., 1987). Sin embargo nuestros 
resultados no se ajustan a esta hipótesis. Estos 
resultados podrían estar relacionados con la his- 
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toria térmica de las especies andinas (Nespolo 
& Rosenmann, 1997) y con la historia evolutiva 
de los taxa relacionada a climas fríos o templa¬ 
dos (Nespolo et al., 1999; Steppan et al., 2007). 
Además, estudios fisiológicos en algunas especies 
andinas (ej. P. xanthopygus) muestran que mayo¬ 
res temperaturas podrían ser una limitante en la 
distribución de la especie, ya que no cuenta con 
la plasticidad en las repuestas fisiológicas y com- 
portamentales para hacer frente a temperaturas 
elevadas, disminuyendo de esta manera su fitness 
(Sassi & Novillo, 2015; Sassi et al., 2015). 

A modo de síntesis, podemos destacar que los 
Andes áridos se caracterizan por un elevado nú¬ 
mero de especies endémicas, las cuales se ajustan 
a un patrón altitudinal positivo, con una fuerte 
relación especies - área. A escala local, se regis¬ 
tra un aumento de riqueza a medida que aumen¬ 
ta la altura y una mayor abundancia en alturas 
intermedias. Los factores causales del patrón de 
diversidad son principalmente el clima y la topo¬ 
grafía, cuyo efecto conjunto explican un 35% de 
la variabilidad de la riqueza. 

La relación inversa entre temperatura y ri¬ 
queza es un hecho de marcado interés, ya que 
los escenarios del cambio climático pronostican 
aumentos de 2-4° C en la temperatura para la 
región andina (Labraga & Villalba, 2009), y esto 
tendría un impacto directo en la biodiversidad de 
la región. 
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